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This paper presents a Localized Radial Basis Functions Collocation Method (LRBFCM) for numerically solving one and 
2-dimensional Fractional Integral Equations (2D-FIEs). The LRBFCM approach decomposes the main problem into 
several local sub-problems of small sizes, effectively reducing the ill-conditioning of the overall problem. By employing 
the collocation approach and utilizing the strong form of the equation, the proposed method achieves efficiency. 
Additionally, the matrix operations only require the inversion of small-sized matrices, further contributing to the 
method's efficiency. To demonstrate the effectiveness of the LRBFCM, the paper provides test problems encompassing 
linear, nonlinear, Volterra, and Fredholm types of Fractional Integral Equations (FIEs). The numerical results showcase 
the efficiency of the proposed method, validating its performance in solving various types of FIEs. 
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1 | Introduction 


COD, iccnsce 


Computational 


In recent years, fractional calculus has gained significant attention for its ability to model complex 
phenomena in various fields such as physics, engineering, and finance [1]—[10], [11]—[13]. Fractional 
Algorithms and Integral Equations (FIEs), a vital component of fractional calculus, have been extensively studied and 
Numerical Dimensions. | applied in numerous applications. To solve these equations, several numerical methods have been 
proposed, including finite difference methods [14], finite element methods [15], and spectral methods 


[16]. Among these methods, meshless methods have emerged as a promising alternative due to their 
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flexibility and adaptability in handling irregular geometries and complex boundary conditions [17]. 
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Radial Basis Functions (RBFs) are a popular choice for meshless methods as they provide a simple 
and efficient way to approximate functions in high-dimensional spaces [18]. RBF-based methods have 
been successfully applied to solve various types of integral equations, including Fredholm and 
Volterra integral equations [19]. However, the application of RBFs to FIEs is still relatively limited. 
Recent attempts have been made to develop RBF-based methods for solving FIEs. For example, 
Dehghan and Mitzaei [20] proposed a meshless method based on RBFs for solving Fractional 
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Fredholm Integral Equations (FFIEs), and Wang et al. [21] developed a similar approach for Fractional 
Volterra Integral Equations (FVIEs). These methods have demonstrated promising results in terms of 
accuracy and computational efficiency. Nevertheless, there is still room for improvement in the 
development of RBF-based methods for FIEs, particularly in the investigation of local RBF methods 
that focus on small subsets of the problem domain [22]. Meshless methods have garnered attention due 
to their straightforward numerical implementation without the need for meshing or re-meshing. These 
methods offer flexibility in handling irregular geometries and can be easily extended to higher- 
dimensional problems. They can employ strong form [23]—[26] or weak form formulations of the 
governing equations [27]—[31]. Local methods have been proposed to address the fully assembled and 
ill-conditioned final matrices obtained by global methods. Popular local methods include the Local 
Multi-Quadric (LMQ) approximation method [32]—[34] and the Finite Collocation (FC) approach [35], 
[36] which reduce the problem into local sub-problems and assemble them into a sparse final global 
matrix. In this paper, we propose a Local Radial Basis Functions (LRBF) method for numerically solving 
PIEs. Our approach combines the advantages of RBFs and local methods, providing a flexible and 
efficient framework for solving a wide range of FIEs. We present a detailed analysis of the proposed 
method, including convergence properties and error estimates, along with numerical examples to 
demonstrate its effectiveness and robustness. The remainder of the paper is organized as follows: 
Section 2 presents basic definitions in FIEs and fractional calculus necessary for the subsequent sections. 
In Section 3, we describe the method of local RBF for discretization of FIEs. Numerical results are 
presented in Section 4. Finally, Section 5 provides a summary of the paper. 


2 | Fractional Integral Equations 


The basic definitions for the Riemann-Liouville Fractional Calculus (RLFC) [1], [37] are presented in 
this section. 


Definition 1 ([37]). Let a finite interval j=[a, b](-0<a<b<)of R. The left and right Riemann-Liouville 
Fractional Integrals (RLFIs) aI,“f and J,“ f of order a € R* are given respectively as below: 


LtO=7 of (t—x)""f(x)dx, (t>a,a>0), (1) 


and 


ItO=Toh (x—t)"'f(x)dx, (t>a,a>0), (2) 


where I(w) is known as gamma function. In the special case of a = n € N+, Eqs. (1) and (2) are the nt 
integrals in the following forms: 


g 1 f n-i 
I 0 rawr (t—x)" Tf (x)dx, B) 


1 
(bl 


Similarly, the below definition for the two dimensional RLFIs of order r is given by considering L'(J) as 


[ox -D £O0d. (4) 


the space of Lebesgue-integrable functions œ: J — Æ" with the following norm: 
a pb 
I wll. = fÈ f I we. yp ll dyax, 
where J = [0, a]X[0, b]. 


Definition 2 ([37]). The left-sided mixed RLFI of order r = (r4, r2) of the function u(x, y) is defined as 
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r = 1 ae i eti _ +4)t-1 
Iu) x,y) = red. x-s) y-t)" u s,t)dtds, (5) 


where u € L! J). r € (0, ©) x 0,00) and © = (0, 0). So, we have 


I. (Su) x,y) = u x,y). 
Il. (Gu) x,y) = N pu x,y)dtds where t = 1,1) and forall x,y) € J. 
Il. (Qu) x,0) = Igu) 0,y)=0 forx €[0,a], y € [0,5]. 
IV. Let’, w€ —1,00) then ( o xy”) = PD Ot) yd yore forall x,y) €J. 


I A4trnjt)PF w+rz+1) 


The interested reader can refer to [38] for more details of the left-sided mixed RLFI. The following is the 
FVIE: 


u(x) = f x) + va Í Kx, t) x-t)*"u t)dt, 0<x<1, (6) 


where a> 0 is a real number, the right-side function f is given, and K(x, t) is the kernel. a =1 is 
corresponding to the ordinary (non-fractional) Volterra equation. Note that when a is non-integer, the 
term (x — t)“ adds only up for x > t. Hence, defining the FFIE as the following form is reasonable: 


u x) =f x) +t} S K (x, t)(x - t)®™tu pat+ f K x, t)( eo) lu tdt]. (7) 


Two-Dimensional (2D) FVIE is defined as follows: 


= 1 < Z — e171 _ 4)\to-1 
u9 =8xy)+ ro I S. x- s) y- t)? K(x y,s,t,u s, t))dsdt. (8) 


3 | Discretization by LRBF 


In this section, the local RBF method is used as a technique for approximation of FIEs. The function ®: 
R+ — R is RBF which is defined as the function of distance r = Ilx — xl [39]. Our choice of RBF is the 


MQ (r) = yr? +c? which belongs to a class of infinitely differentiable global RBFs. Consider the 
following VIE of fractional order: 


u(x) =f(x) + Tob K(x, t)(x —t)*‘u(t)dt, O<a<l. (9) 


A set of N distinct points X= {x,X9,....,Xy} in R, which are called centers is used to discretize the domain 
of the problem. There are not any restrictions on the location of the centers or on the shape of domains. 


The solution u can be approximated at each of N centers by a localized formulation as below: 


ux,)* ix)= $ AS (IK. 2 xsi), (10) 
j=l 


where n is known as the number of nearest neighbouring points {xF} which surrounds the collocation 
point x, , containing the collocation point itself. The stencils are center and its n-1 neighbours. 4; 's are the 


unknown coefficients and @ is an RBF. By applying the interpolation conditions 
n 
u(x?) = ü x8) = X Ajo (Ix: E xii] Se eee (11) 
j=l 


On each N stencil we have an nxn linear system @A5 = us. The matrix @ with elements Pij = (IG - 


xl) is called the interpolation matrix. If (x) is a positive definite RBF and all collocation points are 


distinct, then the interpolation matrix of RBF @ is non-singular and we have 


A8 = 1uS, (12) 


where 

AS = AG, Aj,.-, AR), ue = Uj, u5,..., ug). 

Therefore, the estimated solution a(x, ) can be approximated by the given nodal values u (<) at stencil 
points: 

ü = OA = oh lus = Yus, (13) 
where f = (IIx, - xi). If we rewrite Eq. (13) in the approximate solution terms U(x;) at all collocation 
points, then 

UCx,) = Yu, 
where W is a sparse matrix of order NXN which has at most nxN nonzero elements. Considering Eq. 
(6) and Legendre-Gauss-Lobatto nodes and weights, we have 

a-l 


y= (== Yo, (x,-t,) K(x,, t, u"®(t, ), j=1,2,....N, i=1,2,....n, (14) 
T(a) k=1 


where w, and t for k =1, 2, ..., m are Legendre-Gauss-Lobatto weights and nodes, respectively. Thus, 
we have 


Aŭ =B. (15) 


As a system of linear sparse equations. The approximate solutions at all collocation points can be found 
by solving Eg. (15). In the same manner, for 2D-FIE we have 


up) = X Sy (lp = pyll) =Clo p => Cao, (16) 


where p=(x, y) and p, = (ep yy) € R?. We can remove the dependence on the RBF expansion 
coefficients from Eq. (16) by the following: 


up) = $) coips- pill) = O°C* = b8tus = War. a7) 


Substituting Egs. (16) and (17) in Eq. (9), we have 


1 x gY 
Sa — ey-1 — +)r2-1 -1 ñ a T8 
Wus = 8 xy) +E Í Í xs) y-H21K(x,y,s,t.0'o s, Dü)dsdt. (18) 


Substituting the given collocation points into the above equation and applying Legendre quadrature 
integration formula, we obtain 


1 m m 


tistan) prg LEME) (ya) Keyta) 


4 | Numerical Experiments 


In this section, we apply the proposed method on some test problems and then we present the numerical 
results. For evaluating error estimation, the quantity M = 2J can be defined in which J is the maximal 
level of resolution and the interval [a, b] can be divided into 2M subintervals with the equal length, and 
N = 2M + 1, where N is the number of collocation point xı. When the problem u = uex(x) has a known 
exact solution, the differences 


ADi) —u,,(x,), 1=1,2,...,2M +1, 


can be calculated and the error estimation can be estimated as 


Ou = max, [Ax (1) | 


or 
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Uexll 
2M" 


When the problem has unknown exact solution, the problem can be solved by some level of resolution J, 


Oex = lŭ ~ 


and u,(x) represents the result. Then, these calculations are repeated for J+1 getting uy,,(x). The differences 
A, (xı) =u; (x,) Uy, (x), 


are defined, where xi, 1 = 1, 2, ..., 2M + 1 are the collocation points at the level J. The error estimations is 
then obtained with 


6, =max|A,(x,)|, œ= A EDIM). 


4.1 | Test Problem 1, Abel’s Integral Equation of the Second Kind 


The below Abel’s integral equation of the second kind is given [40]. 
1 2aresin (Vx ) u(s) 
u (x) = + = f 
1+x 1+x o VS+X 


dt, 0<x <1. (19) 


Table 1. The RMS error obtained with different values of 


nodal N for Example 1. 
LRBF GRBF 
N n=6 n=9 n = 16 
20 2.6 x103 2.2 X 105 2.1 x 103 2.2 x 1073 


30 2.1 x 10% 2.2 x 103 2.5 x 103 2.3 x 103 
50 2.1 x 10% 2.4 x 10% 1.2 x 10% 2.2 x 10% 
85 1.9 x 10% 3.0 x 10% 1.7 x 102 1.6 x 10% 


Fig. 1. RMS errors for local and global RBF methods for Example 1. 


Which its exact solution is u(x) = + . Numerical results versus the numbers of nodal points and stencil 
are shown in Table 1. Fig. 1 shows the error curves for local and global RBF methods with 50 number of 
nodal points. Fg 2 presents the RMS error with 50 nodal points and versus stencils. According to Figs. 7 
and 2, the optimal choice of the stencil size for 50 nodal points is n = 6. The local method with n « N is 


often just as accurate as the global method. 


Fig. 2. RMS error for N = 50 and versus stencils for Example 1. 


4.2 | Test Problem 2 


Consider the equation 


u(x) — =|] @ -x —t)(x — t)" u (t)dt + ine -x —-t)(x -t)" u (t)dt]=x° + x 


T(a) 


15 


When « = 1, corresponding the non-fractional equation, u(x) =x? - 3, x €/0,3] is the exact solution. 


Results obtained with N = 45 nodal points and stencil size n = 15 are given in Fs. 3 and 4 and Tabh 2. 
The analysis of error has been done for « = 0.9, Table 2 shows the results. Approximate solutions with 


various « are presented in Figs. 3 and 4. Considering the results, accuracy of the global method is more 
than local method, but the local method is faster. The results are the same as that of [41]. 


Table 2. Numerical results for Example 2 with « = 0.9. 


LRBF GRBF 
J 2M oy OF by Oj 
3 16 0.4784 2.9 x 107 0.0329 2.1 x 107 
4 32 0.2482 7.8 x 10° 0.0325 1.0 x 10-7 
5 64 0.1772 1.4 x 107 0.0367 5.7 x 10% 
6 128 0.8643 6.8 x 10° 0.1239 9.6 x 107 


t, 


Ba 


ao 
w l 
SrerevErri A 


+ 


+ 
=i 
— 


ape 


Fig. 3. Numerical solutions of Example 2 for various «. 
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4.3 | Test Problem 3, 2D Nonlinear FVIE 


Consider the following 2D nonlinear FVIE: 


u(x, y) — —+,| 
n 


7 


T 2 2 r =f 42 
f, f (x -s)? (y-t)? /xytfu(s, t)] aas |= ys ix y? za (20) 


Fig. 4. RMS error versus N for Example 2. 


yoxy 
3 
3. The RMS error versus the number of the nodal points N with 3 x 3 and 5 X 5 stencils is given in Fig. 5. 


Fig. 6 shows the error with N= 82, 3 X 3 and 5 x 5 stencils. Figs. 7 and 8 show the cross section of the 


For which the exact solution is u(x, y) = —. The results with shape parameter c = 4 are presented in Table 


approximate solutions and errors with N = 102, 3 X 3 stencils and fixed values of y. 
4.4 | Test problem 4, 2D nonlinear FFIE 
Consider the following 2D nonlinear FFIE: 


1 T eT 7 1 3 
u(x, y) -l J, T-9)?T-9?54s(y - x Iuls, 0] aa] =f(x, y), 


roe) 
2 2 
where 
T=1, 
ne go ee oles (21) 
322560 l 
Table 1. Numerical results obtained with different number of nodal points for Example 3. 
3 X 3 Stencil 5 x 5 Stencil 
N 5ex RMS Sex RMS 
6? 13x107 5.8480 x 107° 8.5 x 107 3.5273 x 1073 
T 3.0 x10 1.2000 x 1073 4.6 x 1072 2.9653 x 1073 
82 1.1 x 107? 3.7382 x 1074 3.0 x 107° 1.0000 x 1073 
92 1.8 x 105 2.7669 x 1074 2.9 x 103 9.5724 x 10-4 
102 2.1 x 10-4 8.6700 x 1075 5.1 x 10-3 3.3647 x 1074 
112 3.7 x 10-4 9.2053 x 1075 4.2 x 1072 3.2812 x 10-3 
122 7.4 x 1073 8.1000 x 1074 7.5 X 107? 1.4528 x 107? 


Table 4. Numerical results obtained with different number of nodal points for Example 4. 


Fig. 5. RMS error versus N for Example 3 obtained with 
3x3 and 5X5 sencils. 


Fig. 6. RMS error obtained for Example 3. 


= 


3 x3 Stencil 


5 xX 5 Stencil 


N ĉex RMS ĉex RMS 

6&2 12x103 5.9600 x 10% 8.2 x 103 8.9394 x 104 
7 10x 10% 5.8873 x 1074 6.6 x 103 5.6580 x 10% 
s2 1.0 x 10% 2.3020 x 104 5.4 x 103 3.4530 x 104 
92 664x104 2.9741 x 105 3.2 x 103 3.0230 x 104 
102 3.4 x 103 1.0102 x 104 5.8 x 104 2.9436 x 10° 
112 4.3 x 103 5.6420 x 10°4 6.3 X 104 2.3549 x 104 
122 6.3 x 102 2.3751 x 103 2.6 X 10°2 1.1205 x 1053 


nodal points are shown in Table 4. 


. . a ss 1 . . 
And the exact solution of this equation is u(x, y)= x? Be y. The results with different number of 
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Fig. 7. Cross section of approximate solution obtained with N = 
10? and 3X3 stencils for Example 3. 


4.5 | Test Problem 5 


Consider the 2D nonlinear FVIE: 


u(x, y) — ronal -s)3(y- p(s, hey [u(s, t)f aas =f(x,y), 
3 3 


0 
where 


19683x° (77y +18y° (3y -7)) 
22422400 l 


f(x, y)=x’y y-l 


B(., .) denotes the two dimensional Beta function. Exact solution of this problem is u(x, y) = x?(y 2-y). 
Table 5 illustrates the numerical results with different number of nodal points for this test problem. The 
results obtained in Test Probkms (3), (4) and (5) with local RBF method are more accurate than the results 
obtained in [37]. 


[tt 


Fig. 8. Cross section of approximate solution obtained with 
N = 10? and 3x3 stencils for Example 3. 


Table 5. Numerical results obtained with different number of nodal points for Example 5. 


3 X 3 Stencil 5 X 5 Stencil 
N Sex RMS Sex RMS 
2 5.8 x 1073 8.9341 x 107-3 1.7 x 10-2 1.2324 x 10-3 
2 5.7 x 1073 1.9341 x 10-3 1.0 x 10-2 1.0561 x 10-3 
2 5.3 x 10-3 7.9201 x 10-4 8.9 x 1073 3.7213 x 10-3 
2 8.3 x 10-4 3.2047 x 1074 6.3 xX 10-3 7.7001 x 10-4 
102 5.8 x 1074 2.6000 x 10-4 2.3 x 1073 9.1540 x 10-4 
112 3.2 x 1073 3.9805 x 1074 9.3 x 1074 1.5170 x 1074 
127 2.5 x 10-2 1.2005 x 1073 5.3 x 10-3 5.3072 x 10-3 


5 | Conclusion 


A meshless local RBF method was proposed for solving one and two-dimensional FIEs. In the local 
RBF method, the approximate solution at stencil centers was expressed in terms of the given nodal 
values u(xj), which correspond to the n-nearest neighboring points. This allowed for the determination 
of the approximate function values at the nodal points without the need to calculate unknown 
coefficients Aj. Numerical results demonstrated that the local RBF method outperformed the global RBF 
method, particularly in two-dimensional FIEs. This suggests that the local RBF method is more suitable 
for high-dimensional problems. 


The effectiveness of the proposed method can be attributed to the following reasons: 


I. The use of the strong form equation and collocation approach simplified the method. 
II. By employing the localization approach, the matrix operations only required the inversion of small- 
sized matrices, resulting in a sparse final global matrix. 


Based on these advantages, it is recommended to utilize the proposed method for tackling more complex 


and similar applied problems. 
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